POSITIVE TRIGONOMETRIC QUADRATURE FORMULAS 
AND QUADRATURE ON THE UNIT CIRCLE 



FRANZ PEHERSTORFER*^ 



Abstract. We give several descriptions of positive quadrature formulas 
which are exact for trigonometric - , respectively, Laurent polynomials 
of degree less or equal n — 1 — m, < m < n — 1. A complete and simple 
description is obtained with the help of orthogonal polynomials on the 
unit circle. In particular it is shown that the nodes polynomial can be 
generated by a simple recurrence relation. As a byproduct interlacing 
properties of zeros of para-orthogonal polynomials are obtained. Finally, 
asymptotics for the quadrature weights are presented. 



1. Introduction 

I 

Let / G No := N U {0}, 7 G {0, 1/2}, and let 7i,^ = { E Ofc cos(fc + 7)99 + 

fc=0 

bk sm{k + 7)(/9 : ak, bk G M}. By 71 := TI.o we denote the set of trigonometric 
polynomials of degree less or equal /. Further let a be a positive measure on 
[0,27r] normalized by Jq^ da = 1. We call a quadrature formula (qf) of 
the form 
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r2TT " 
(1) / t{ip)dai^) = fisti^s) + Rn{t) 



with < Lfi < < ■■■ < fn < 2vr, ni, ...,fin G M, a trigonometric (n — 1 — 
m, n, da) qf if i?„(t) = for t G Tn-i-m] n — 1 — m is the so-called degree of 
exactness. Note that for n G N, n = 2(n-|-7), 7 G {0, 1/2}, the trigonometric 

n 

—i{nip+ J2 Vs)/'2 

nodes polynomial 11^=1 sin(^^^) = e ''=1 nr=i(e*^ " e*'^»)/2i is 

from Tn,-y 

Obviously qf ([l]) is a (n — 1 — m, n, da) qf if and only if 
(2) ^ye-^^^daM = 5^;U,e-^^^= k = 0, . . . ,±{n - 1 - m) 







s=l 
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which is a qf on the unit circle exact for Laurent polynomials dkz'^. 
d/c S C of degree < n — m — 1, which we call (n — 1 — m, n, da) Laurent 
qf also. In other words we are looking for a characterization of quadrature 
formulas which integrate exactly the first n—l — m moments ci , . . . , Cn-i-m 
of the Caratheodory function 

(3) F{z) = i- / "^M^) = 1 + 2 E 

^==1 

Because of convergence and stability reasons it is of upmost interest to 
have positive (n — 1 — m,n,da) qf, i.e., that all quadrature weights fig are 
positive. Such qf will be studied in this paper. So far it is known that qf with 
highest possible degree of exactness, that is, (n — l,n,da) qf have positive 
weights. These (n — 1, n, da) qf can be considered as the counterparts of the 
Gauss qf for algebraic polynomials and are called Szego qf nowadays. The 
Szego qf can be described as follows (see [HI pp. 14-16], in particular (10.10), 
but also [H [131 [221 [2S]): fi,---, V'n, < <...<(/?„ < 2tt, generate a 
(n — 1, n, da) Laurent qf if and only if there is a r/ G T := {z G C : |2:| = 1} 
such that 



(4) Z^n-l{z)+vK-liz) = ll( 



n 

' Z 

s=l 



Here and in what follows ^n{z) = + . . . always denotes the monic poly- 
nomial of degree n orthogonal with respect to the normalized measure da, 
i.e., 

27r 

(5) J e-'^'^^nie''^)da{ip) = fc = 0,...,n-l 



It is well known that the 'I'ri's satisfy a recurrence relation of the form 

(6) ^n{z) = Z<!>n-l{z) - an-lK-l{z), 

with \aj\ < 1 for j = 0, 1, 2, . . . and where $*_i(2;) = z"~-'^¥„_i(i); hence 

(7) $:(z) = K~liz) - an-lZ'Pn-liz), 

Furthermore let us recall that for any k E T 



(8) Z^n-l+K^*n-l = ll('-e'^^) 

s=l 

where < ipi < ■ ■ ■ < ipn < 27r, , since ^n-i has all zeros in B = {\z\ < 
1}. Following ^13j a polynomial of the form ([8|) is called para-orthogonal 
polynomial. 
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The Szego quadrature weights are given by, Zg 



where ^„ is the so-called polynomial of second kind, i.e. 

(10) ^n{z) = -j J3^(^n(^) - ^n{e^n)d<y{^) 



which satisfies the recurrence relation 

(11) ^'n(^) = Z^n-l{z) + an-lK-l{z). 

Note that a„_i in ([U]) is just replaced by — a„_i. Recently Szego quadrature 
formulas have been studied and applied by several authors [31 [6l [TTl [TOl [H]. 
For para-orthogonal polynomials see O [Ml HZ] ■ 

As mentioned already, here we study positive (n — 1 — m, n, da) qf. There 
are two possible approaches, one approach comparable to the algebraic poly- 
nomial approach and the other one via orthogonal polynomials on the unit 
circle (abbreviated by OPUC). Each of which has its advantages as we shall 
see. In the next section we derive some characterizations without using 
OPUC's. In fact what concerns Theorem 12.21 we even do not know how to 
obtain it via OPUC's and what concerns the other properties as the orthog- 
onality condition given in Lemma 12.11 they follow in a more natural way 
than with the help of OPUC's, compare Lemma 13.51 and Corollarv 13.61 In 
section 3 we give a simple complete description of positive trigonometric - 
resp. Laurent qf with the help of OPUC's. In this problem the powerfullness 
of the approach via the unit circle turns out nicely. We note that section 3 
and section 4 can be read without having studied section 2. 

2. Basic characterization 

First let us take a look at the degree of exactness. Similarly as in the 
case of algebraic polynomials the degree of exactness is connected with the 
orthogonality property of the nodes polynomial. For Szego qf see [3] and 
for the associated orthogonal trigonometric polynomials and their five term 
recurrence relation see jl2j . 

Lemma 2.1. Let n G N, n = 2(n + 7), 7 G {0,1/2}. A qf of the form 
([1]) is exact for t G Tn-\-2y+h ^ £ No, if and only if the trigonometric nodes 
polynomial Tn^^{ip) satisfies the orthogonality conditions 

(12) / tiip)Tn,^i^)da{ip) = for t G 71,^ 
Jo 

Proof. Case 1) n = 2n, i.e., 7 = 0. 

Sufficiency. First let us consider the case I = 0. Let ^ G [0, 2tt) be no 
zero of Tfi and arbitrary otherwise. Since the Lagrange trigonometric poly- 
nomials with respect to the nodes fi, V2n,C a-re polynomials from Tn the 
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interpolation quadrature formula 

(13) / t{(p)da = y2 Hst{(Ps) + fi2h+it{0 

holds for every t G Th and is unique. Now any q G Ta can be represented in 
the form 

(14) qi^)=cTni^)+ti^) 

where t G 7s is the unique polynomial satisfying the interpolation conditions 
= t{ips) = qi^s) for s = l,...,2n and c = q{(,)/Tn{S^). Hence 
applying (fT3]) to q, taking into consideration the orthogonality property of 
Tfj, we obtain 

^ f-2iT 2 /■2t I f2n 2n 2n 

IT <lda= — cTnda + — / tda = V/ist(9?s) = ^fisqi^s) 
27r Jq 27r 27r ^ ^ 

In case of a higher orthogonality property of T^jwe represent g G Tn+i 
the form 

g = vTn + t 

where t is as above and v £ Ti and the assertion follows as above. 

The necessity part follows immediately by applying the quadrature for- 
mula to tTfi. 

Case 2) n = 2n + 1, i.e., 7 = 1/2. 

Sufficiency. Any q S Tf^+i+j can be represented in the form 

Q = ^'«,l/2^ri,l/2 +* 

where t £ Th interpolates q at the 2n + 1 zeros of and where f;,i/2 £ 

T'i,i/2- Applying the unique (n, 2n + l,da) quadrature formula to q, 
putting ^ = V32n+ij the assertion follows by the orthogonality property. 
Necessity is immediate again. □ 

We note if / < n — 1 then the orthogonality property (|12p does not imply 
that Tfi^^ has all zeros simple and in [0, 2tt) there may appear also double or 
complex zeros. 

To find a criterion on the nodes polynomial such that the quadrature 
weights are positive is more involved (in contrast to the algebraic polynomial 
case, where the corresponding counterpart can be derived in one line, see 
e.g. m (1-3)] 

Theorem 2.2. Let n e N \ {1}, n = 2(n + 7), 7 G {0, 1/2}, and let the qf 

1^ be a (n + 27, 2(n + 7), do") qf. Denote by T{ip) the trigonometric nodes 
polynomial, which is from Tn,^, and put 

(15) 5(V) = ^ r cot (^^) (T(V) - T{ip))da{ip). 
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Then S{ijj) E Th,'y and 

(16) Us = — T-n — r for s = 1, n. 

In particular, all quadrature weights fig are positive if and only if S{(p) has 
n simple zeros 6j, < 9i < 02 < ... < On < 27r, which separate the zeros 
ifj ofT{ip) such that ipi < 9i < ip2 < ■■■ if sgn{ST){0~^) > 0, respectively, 
ei<ipi<e2< ... if sgn{ST){0+) < 0. 
Furthermore, setting 

Tn{e'n = e^t'^r(^) and L^„(e*'^) = ie*t^5((/p), 

the ratio of the polynomials ijJn{z) and Tn{z) has an expansion of the form 

where Ck = J^^ e~^^'^da{(f) and h + 2'y + 1 is the precise degree of exactness 
of the qf dH). 



Proof. Let R{ip) G Th,-y be the polynomial such that, z = e*"^, 

ie'^'^Rjip) _ -A z+_e^ 



(18) ^S^ = Em. 



Putting Unie'^P) = ie'2'PR{ip) and T„(e*'^) = e^2'^r(99) we may rewrite 
equation (jlSp in the form 

(19) ^-i = -2y 

^ ' s=l 

where we used the fact that 



s=l '' 



since the qf is exact for constants and a is normalized. Note that i^n(O) = 
— T„(0) or in other words that 

(20) 

n n n n 

T {if) = On cos —ip + bnS\n—Lp + ... and i?(99) = — 6„ cos — (/9 + sin — 99 + ... 

Furthermore it follows from (|19p that 



lips 



2/ise"^« 



which yields that 

(21) ^-^(y^^) = e-^t^-u;„(e^^O ^ Unje'^^) ^ 2^ 
§i^s) ^{e-^^^'rn{e^^))^=^^ ^(^^(e^^))^^^^ i 
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The goal is to show that 
(22) S = R. 



Let us prove ()22p when n is even. 
Using the facts that 



cot ( ^ ) (cos ^ — cos Lp) = sin + sin ip 



and that 

, , cos kih — cos kip 

(23 ^7 = 2 cos /c + ... 

cos ^ — cos 

is a cosine polynomial of degree k — \ va ijj with leading coefficient 2, it 
follows by straightforward calculation that for /c G N 

cot I ^ ) (cos /cV ~ cos A;(/3) da{(p) 
V 2 y 

(24) = / cot ( ^ — ^ ) (cos V' — cos w) f —] da(ip) 

Jo V 2 y \ cosil^ - cosp J 

2tt 



da) sin kip + t('0), 
where t G Tfc-i- Furthermore, 







cot ( — ^ — j (sin kip — sin kip)da{ip) 

■ I f^"" f ip - f sin kip sink(p\ 
smip / cot — - — — — ■ : da{p) 



(25) 







2 J \ sin^/^ sinc^ 



f'^^ / — ^ \ . . , . , sin few , , , 

+ / cot (sm — sm (/?) da[ip) 

Jo V 2 y sm^J 



2tt 



= — ( / do") cos /cV" + i(V') 

JO 

where t G 7fc-i, and where in the last equality we used the facts that 
sinfex/sinx = 2 cos (A; — l)x + ... is a cosine polynomial of degree k — 1 
in X and thus (j24p can be applied to the first integral and that the second 
integral is from 71. Thus S{'ip) G Tn and by ()20p it follows that S and R have 
the same leading coefficients, i. e., i? — 5 G Th-i- Writing T{ip) in terms of 
sin {■^-^^) it follows that the integrand from ()15p at = tps is from 7fi and 
thus we may apply the quadrature formula to the integrand and obtain 

(26) S{ps) = -^J's'^T' (ps) = R{^s), s = l,...,2n, 



where the last equality follows by (j2ip . Hence S = R and the assertion is 
proved for n even. 
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For odd n we consider 

^^"^^ = ii''^'"* (^) ^^'^^^^'^^ " {qTKip))da{^) 
where q is any polynomial from 7o^i/2 with q{ip) > on [0, 27r). Since qT £ 
Tfi+i it follows by (p^ and ([25]) that S G 7fi+i. Furthermore, we may apply 
the quadrature formula and obtain that 

(27) S{ips) = -fis2q{iPs)T' (ifs) = q{ips)R{^Ps) for s = 1, . . . , 2n + 1. 
Next we claim that 

(28) S{^P) = qWSii^) 
To prove ()28p it suffices to show that 

cot (^y^) {qW - q{f)) n^)da{^) = 

which follows immediately by g(V') — <z(v) = const sin (^^y^^ and the or- 
thogonality property of T with respect to 7o,i/2- 

Therefore, by (|28p and (|27|) . if we are able to show that the leading 
coefficients of S{ip) and {qR){ip) coincide, i. e., S — qR € Tn representation 
([T5|) follows for n odd also. 

Because of (pOj) we know the leading coefficients of {qT){ip) and {qR){ip). 
By applying ([M]) and ([25]) to iqT){'ip) - {qT){ip) it follows that the leading 
coefficients of S are equal to that ones of qR. Hence S = qR and the theorem 
is proved. 

The statement on the positivity of the fig follows immediately by (126p . 
The last equality in (jl7p follows by partial fraction expansion and ([2]). 

□ 

Theorem 12.21 (which turns out to be crucial in Szego-Kronrod quadrature 
|21j for instance) does not show how to obtain a constructive description of 
positive trigonometric qf. Such a description will be obtained by OPUC's 
in the next section. 

3. The Characterization Theorem 

First let us recall that the monic polynomials ^niz) of degree n, n G No, 
are the polynomials orthogonal with respect to da, where ^ Jq^ da = 1, 
and that the Oj's are their recurrence coefficients, see 

Theorem 3.1. The following four statements are equivalent: 
a) 

(29) [ e-'^'^dai^) = V ^u.e-*'^^^ for = 0, . . . , ±(n - m - 1) 

^=1 

where <(/?!< (^92 <•••< V'n < 27r and fj,i, . . . , fin G M+. 
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b) There exists a polynomial ^n-i{z) = -z" + ... generated by the 
recurrence relation 

(30) ^j{z) = zl>,_i(z) - a,_il>*_i(z), 
where 

(31) |aj I < 1 for j = 0, . . . , n — 2 and aj = aj for j = 0, . . . , n — m — 2, 

and a r] gT such that 

n 

(32) z^n-i+vK-i = ll^^-e'''') 

s=l 

c) There exists a polynomial qm{z) = z"^ + • • • which has all its zeros 
in B and a r] gT such that 

n 

(33) Zqm^n-n,-l + = " ^*'^') 

s=l 

d) There exists a polynomial pn-i{z) = z"'^'^ + . . . which has all its zeros 
in B and a r] gT such that 

n 

(34) zpn-i+VP*n^i = lli^-e''^^) 

s=l 

and 

if: n— m— 1 

(35) ^^r' " =1 + 2 E c,/ + 0(."-'") 

27r 

where Ck = J e~^'"^da{(p) 



Proof. a)=^ b) Since ^i, . . . , ^„ > it follows that the sequence 

n 

(36) Cfe = E ^.e-*^'^= A: = 0, . . . , n - 1 

s=l 

is positive definite, see [I]. Thus the polynomials j = 0, . . . , n — 1, 
orthogonal with respect to the positive definite sequence cq, . . . , Cn-i satisfy 
a recurrence relation of the form (I30p with |aj| < 1 for j = 0, ... ,n — 2. 
Since by assumption 

(37) ck = — [ e-'^'^dai^) k = 0, . . . ,n - m - 1 

2vr J 

it follows that, see e.g. j8, (3.2)] 

(38) dk = CLk for A; = 0, . . . , n — m — 2 
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Finally we know [8l Thm.6.1] that for any a^-i with |a„_i| = 1 

- ^^"-^ " = 1 + 2 V c./ + Oizn 

z$„_i +a„_i$;_i ^ 

(39) „ . 

s=l 

Putting a„_i = -e*'^i$„_i(e*'^i)/'l*_i(e*'^i), that is, the denominator in 
(pl) vanishes at e*'^S it follows that, see [8, (9.5)] 



(40) zl>„_i +a„_il>;_i = ]J(z-e^ 



b) =^ a) By (pOj) and the first condition of (f3T]) 4>,„_i(z) is orthogonal with 
respect to some positive measure a. Thus (j32p holds for any r/ with |?7| = 1 
and 

n 

(41) — /" e-^'='^da(vp) = V ^ise-'^'^' for /t = 0, . . . , n - 1 

— 77 

where /ii , . . . , /i„ G M"*" . Since = Oj for j = 0, . . . , n — m — 2 it follows that 
$j = $j and ^'j = for j = 0, . . . , n — m — 1 and thus the first n — m—1 
moments of a and a coincide, that is, 

(42) / e-''"^dd{^) = [ e-''"^da{<f) k = . . . ,n - m - I 



which gives by (I4ip the assertion. 

b) =^ c) Let 

(43) qj{z) = zqj_i{z) - r]an-i-jq*j^i{z) j = 1,2,... ,m 

with qQ{z) = 1 then it follows by induction arguments using ([6]) and ([7]) that 

(44) z^n-i + ri^*n-i = zqj^n-j + VqjK-j 

In view of (I43p qj has all zeros in ID, see e.g. [22]. 

We note that it follows analogously by induction arguments that 

(45) z^n-i - r]K~i = zqj^n-j - VQjK-j 

c) =^ b) Since qm{z) has all zeros in D q^ can be generated by a recurrence 
relation of the form 

(46) qj{z) = zqj^i{z) - rjbj^iqj^iiz) j = l,...,m 

where \bj\ < 1 for j = 0, ...,m — 1 and qo{z) = 1. Hence defining the 
polynomial ^n-i by the recurrence relation ^ with recurrence coefficients 
for OPUC's ao, . . . , a„_m_i, 6m- 1, . . . , 6o it follows as above that (f32]l holds 
which implies b). 
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c) d) We put 

^^^^ Zqm'^n-m^l + r]q*mK-m-l = ^Pn^l + VP*n~l 

that is, 

2zpn-l = zqm{^n-m-l + 'i'n-m-l) + mm (*^n— m— 1 ^n— m— l) 

(48) 

Let us demonstrate that Pn~i has no zero in D. It is known, see (1.3.82), 
(3.2.62) and (3.2.63) from t22j that + 'i'n^m-i has no zeros on B. 

We claim that 

(49) z{<Pn~m^l-'fn-m-l) .^ for Izl < 1. 

n—m—l I n—m—1 

By the maximum principle it suffices to show ()49p for \z\ = 1. Since on 
\z\ = 1 

I z{^n-m-l ~ ^n-m-l) i _ i ^n-m—1 ~ ^n-m—1 i 

V ^ I 1 + ' ~ ' $* 1 + * * 1 

n—m—1 ' n—m—1 n—m—1 ' n—m—1 

and since by [H (10.3)] 

(51) Re ^r'"~ ■ > 



it follows that the last expression in (jSOp is less than 1 and thus (09]) is 
proved. S ince qm has all zeros in O there holds |(^^(z)/^^(2:)| ^ 1 on ID), 
which implies by (j49p and (j48p that has no zero in B. Finally, relation 
(j35]) follows with the help of the relation, see [22, Prop. 3.2.9] 

(52) ^n-m-l{z)F{z) + ^n-m-l{z) = ©(z"-"*-^) 

and 

(53) K_^_,iz)Fiz) - K-m^iiz) = Oiz--n 
and dTTl) . 

d) =^ c) The proof runs as in [171 P- 939]. To make the paper selfcontained 
we sketch it. According to O Theorem 18.2]) there exists a function /i : C — )• 
C which is analytic in B and satisfies the inequality \h{z)\ < 1 for z G B, 
such that 

VPn-li^) - ZPn-liz) Z^n-l-miz)h{z) - ^l_^_^{z) 



(54) 



VPn-liz) + ZPn-liz) Z^n — 1—m 



m— 1 



Z 



Isolating h from this equality we find that h can be represented by 

(z) ' 
(55) /i(z) = Mfl, whereg/(z) = "rT(z-z,), Zi£B, 

qi has real coefficients, and I > m. 
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Let us assume that / > ni. Then it follows froni (I54p that — n—i—m 

h + 

"^n-m-i ^iid z^n-i-mh + <I>*_m_i have (/ — ui) common zeros on \z\ = 1, 
which implies that 

(56) + 

n— m— 1 n— m— 1 

has (/ — m) zeros on |z| = 1. But this is impossible, since (see [8, p. 4]) 

(57) nn~m-lK-m-l + ^n-m-lK-m-1 = Kz^'^'^, 

where G R+. 
a) d) Let 

(58) ck = — [ e-'^'Pdai^) = V ^l,e-'^'P' for = 0, . . . , n - m - L 

'^-^ 
Hence, by expansion at z = 0, 

n , n—l—m 

(59) -E^«737I^ = l + 2 E c.z'= + 0(."-'")- 



s=\ k=l 



Now the left hand side is a rational function with real part zero on \z\ = 1 
whose numerator and denominator is a polynomial of degree n. Thus it may 
be written as 



(60) ^Pn(^) -Pnjz) ^ _ y- ^ + e"^^ ^ ^ 

Wniz) + Pn{z) fr{^'z-e'Vs 

hence 

(61) Pn{z) = zpn-l{z) 

Moreover by ([60]) . partial fraction expansion and the suppositions on /i^, 
using also z-^f{z) = — ^^/(e*'^), we obtain that 

(62) = (',P;,-Pn)(e'--) ^ „ 

Thus the trigonometric polynomials lm{r]^^^'^e~^^'^pn{e^'^)} and 
Re{ry~^/2e~*2'^p,i(e*''')} have strictly interlacing zeros, which implies by the 
argument principle that Pn has all zeros in D. By (j60p and (|59p the impli- 
cation is proved. 

d) =^ a) Since := zpn-i has all zeros in D it follows that ?7p* ±Pn has 
all zeros on and that the zeros of rjPn—Pn and r/p* +Pn strictly interlace. 
Hence (j62p and thus (j60p holds with > for s = 1, . . . , n. Expanding the 
second expression from ()60p in a series at z = part a) follows. □ 

Obviously Theorem 13.11 could be proved in a shorter way, that is, by 
showing the implications a)=^b)=^ . . . =^>d)=^>a) only. But by proving all 
equivalences separately we get additional insights, for instance: 
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Remark 3.2. In the proof of Theorem 13.11 c)=^b) we have shown that ^n-i 
is generated by the recurrence coefficients ao, ■ ■ ■ , a„_m-i, &m-i, ■ ■ ■ ,bo when 
Qm is generated by 

(63) qj{z) = zqj^i{z) - ribj^iq*_^{z) j = l,...,m 

qo{z) = 1, with \bj\ < 1, j = 0, . . . , m — 1. Note that every polynomial q^ 
which has all zeros in D can be generated by a recurrence relation of the 
form dll), see Thm. 9.2]. 

The ideas of proof for Theorem 13.1b ) and d) go back essentially to the 
author's papers [16l|T7j where positive qf on [—1, 1] have been characterized 
by transforming the problem to the unit circle. For the connection with qf 
on [—1, 1] see Section 5 below. 

Finally we point out that recurrence relation ([30]) . under condition (f3T|) . 
offers a simple, constructive way to generate positive trigonometric (n — 1 — 
m, n, da) qf. 

With the help of Theorem 13.11 we get easily some information on the 
location of the nodes of positive qf on the unit circle. 

Corollary 3.3. Let /, m e No, m < and let z^n^i-i+^^l^i-i = Ilu=i~''(.^~ 
e*'^"), K € T, be a Szego polynomial. Suppose that Tn{z) = ns=i(-^ ~ 6*"^=) 
generates a positive (n — 1 — m,n, da) qf. Then in every interval (-0,^, V'jy+i), 
u = 0, . . . ,1, there is at least one ifg. 

Furthermore, every common zero of Tn and z^^-i-i + i^^*n-i-i ^ 2;ero 
of Kqi{z) — r]q^{z), where qi and r] is given by the representation Tn = 

Zqi'^n-l-l+mlK-l-l- 

Proof. By Theorem 13.11 and the obvious fact that every (n — 1 — m,n,da) 
qf is a (n — 1 — /, n, da) qf for I > m, Tn has a representation of the form 

(64) Tn = zqi^n-i-l + m*iK~i~l 

where qi has all zeros in B. Hence Tnie^^") = if and only if for some E Z 
(65) 

argr/ + (2fc + l)vr = arg ^|§^(e^^0 = arg ^(e^^O + arg ^^(e^^O 

and the condition for a zero of the Szego-polynomial looks similarily. Now 
recall that the argument of a Blaschke product is strictly increasing on 
1^1 = 1 with respect to (p. Thus on arg increases by 27r 

n — l~l 

and, taking a look at the RHS of (1651) . arg increases by more than 

2tt on ['i/'iy, V'i^+i], that is, T„ has at least one zero in {^p^,tp|y^l) and the 
assertion is proved. □ 

For the special case that T„ = z^n-i + V^n-i property given in 

the above Corollary was proved by Simon [2^ Theorem 2.3] by a completely 
different approach. For related statements in the real case see [2} [T8 | [T9 l [25] . 
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For the convenience of the reader let us reformulate parts of Theorem 13.11 
in terms of trigonometric polynomials. 

Corollary 3.4. < (/?! < 952,... <(/?„< 27r, are the nodes of a 

positive trigonometric (n — 1 — m,n, da) quadrature formula if and only if 
there exists a polynomial <l>„_i(z) satisfying (|30p and (I3ip and a G T such 
that, z = e''^, 

n 

(66) r/-^/2^-"/2 Y[{z - e'^^) = 2i?e{z-tr,V2|,*_^(^)} 

s=l 

Proof. The statement follows immediately by Theorem 13.11 a) and b) and 
the obvious relations, z = e*"^, 



rj ^I'^z 2 



(67) 

= 2i?e{z-t,//2c|,;_^(^)}}. 



z~^n-x{z) + ??l>:_i(2)l = 2i?e{z-t+ir/-i/24.„_i(z)} 



□ 



We need the following lemma which can be extracted easily from [26] or 
P pp. 144-145]. 

Lemma 3.5. Let fj, be a positive measure on [0, 27r], /ei n G N, 7 G {0, 1/2} 
and n = 2(n + 7). Let Pn-i{z) = z^~^ + ... be such that 

f-2n 

/ e-''"^Pn-i{e''^)dfi{ip) =0 for k = 0,...,n-2 
Jo 

Then for any r] gT 

(69) / t{^)Re{r]^l^e-''^'PP:_^{e''P)}d^i{ip) = QioTt(^Tn-i,^ 



Jo 

Proof. Relation (j69p is equivalent to 

r.27r 

g-i(f -lT(fc+7))</'p^_^ (e*^)ti;u(v9) 







(70) 

+r/ / e-^(tT('^+T))^P*_i(e*^)d//((^) = Ofor A: = 0,...,n-1. 





Now the first integral is zero by (j68p and thus the second integral also by 
taking the complex conjugate of it. □ 

With the help of the above Lemma we get an explicit weight function de- 
pending on n with respect to which the trigonometric nodes polynomial has 
maximal orthogonality. Furthermore we obtain the orthogonality property 
(jl2p of the nodes polynomial (j66p by an approach via OPUC's. 



Corollary 3.6. Let n = 2(n -|- 7), where n G N and 7 G {0,1/2}. // 
ifi,...,ipn, < ifi < ip2,--- < ^fn < 27r, are the nodes of a positive trigono- 
metric (n — 1 — m,n,da) qf then there exists a polynomial Qmiz) which has 
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all zeros in \z\ < 1 and a 7] £ T such that Re{r]^^^z~2 (^q'^^'^_-^_^)(^z)}, 
z = e*"^, is the nodes polynomial and 

'\{y.)ReW/'e-'^'^^{ql,K-i-Jie'nh 



(71) 70 " "-^-"^^^ " \{q*mK-i-J{e'^)^' 

= for t G 7^-1,7. 

Moreover, 
(72) 



f\{^)Re{r,^'^e-'''^^{ql,K_^_^J{e'^)}da{^) = for t G T^-i-m,^ 
Jo 



Proof. By Theorem \3.1\ q^^„-T-^ has all zeros in D, hence, see e.g. [3 [25], 

Applying Lemma 13.51 the orthogonality property ()7ip follows. 

Since is from 7^ we obtain from (j7T]) . by putting t((/?) = 

^('/^)km(^*'^)P) ^ ^ Th-i-m,y, that ([72]) holds with respect to the measure 
(i(^/|$;_i_„(e*'^)|2 instead of da. But, see [H Thm 2.2, p. 198], 

r ^""^"<'" - r for . ^ 0, 1 - ™ 

and thus relation ()72p follows, taking into consideration the fact that the 
orthogonality property ()72p depends on the first n — m — 1 = 2h + 2-y — m — l 
trigonometric moments of da only. 

□ 

We note that in general Im{z^tryi/2|)*_^(2;)}, given by ([30|) and 

(j3T]) . does not coincide with Imjz^ 2"7^i/2(gi^$*__j^_^)(2;)}, since 

where qm is given by qj{z) = zqj-i{z) + r?a„_i_jg*_^(z), j = 1,2, ...,m with 
qo{z) = 1, which should be compared with (03]). 

4. ASYMPTOTICS OF WEIGHTS 

As a consequence of Theorem 13.11 we obtain the following representation 
of the quadrature weights needed in what follows. Note that only the $n's 
and ^rra's and no ^'„'s appear in the representation. 

Corollary 4.1. The quadrature weights fig from i29\) can be represented in 

the form 

(73) 

^I^n—m~l^s |'?m(2^s)| 

r]^l_^„iq*m){Zs){z^n-rn-iqrn + ri^l_„,_iq*m)' {Zs) 

where = e'^- and K^-m-i = '^U]Z^~\^ " lojf ) 
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Proof. By (P, (jS]) and (US]) it follows that 

2Zsiz^n-m-iqm + VK-m-lQ*myiZs) 

Now, at the zeros of z^n-m-iQm + "n^n-m-ilm have 
where in the last equality we used the known fact that [8, (5.6)] 

(75) ^n-m-lK-m-l + n-m-lK-m-l = i^n-m-l^"""""^ 

which gives by (j74p the assertion. □ 

Theorem 4.2. Let da{ip) = f{ip)dip be positive and from Lipj, < 7 < 1 on 
[a, 13] C [0,27r]. Suppose that ipi^n, ■■■,^n,n, < (^i,n < ^2,n < ■■■ < ^n,n < 27r 
generates a sequence of positive {n—m{n)—l, n, w) qj with quadrature weights 
fj.s,n- Furthermore let us assume that the associated qm{n) satisfy uniformly 
on [a, /3] 

hm f 1 - 2 ^eie^^^kl^^^}^ = g{^) 

Then uniformly for ips^n £ [a + e, f3 — e],£ > 0, 

(76) = -jj r + o(l) 



Proof. First let us note that (I73p can be written in the form 

(77) 
-4^s = 

i^n-^-i|gn^(e^^OP 

Im{e-*t'^»(ryi/2^;;,ci,;_^^_J(e^^.)}^Re{e-*t'^(r/i/2g;;ci>;_^_,)(e*^)}^=^, 

Now 
(78) 

Im{e-t^(r?V2g*^$;_^__J(e^^)}^Re{e-tn^'/'C^;-,n-i)(e^^)} 

-Re{e-t^(ryV2g;^ci,;_^_^)(e'^)}Aim{e-t^(r/V2g;;ci>;_^_^)(e^'^)} 

d 



Im{(e-t^(r?V2g*^$;_„_i)(e^^)) — (e-t^(r?i/2g;;$;_^_,)(e^^))} 
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and thus, since the e*"^" are the zeros of Re{e~*t'^(ryV2g,*^<|,*_^__^)(gj</')|^ -^g 
may replace the denominator in ()77p by the last expression from (j78p which 
yields 



(79) ^isM^K-m-lW^')? 



4imlZ2 ^e (g;;^$;__J(e>^.) I 



Now by the assumptions on / it is known, see e.g. [25^ [22] that uniformly 

on [a + e, /3 — e], e > 0, 

taking into consideration that 

(81) limi^,_^(„)_i/2 = 1/15(0,^7), 
where D{z,a) is the so-called Szego function, that is, 

(82) D{z) = exp{^ I ^ log md^} 



Next let us prove that ([50]) and (jSTj) imply that uniformly on [a + £, /3 — 

e],e > 0, 

(83) -1 : / ■ X 

n'$^„_i_^(e*'^)' n^oo 
Indeed by the local version of Bernstein's inequality 



(84) ((2n - m - 1) max^|cl>*„_„_i(e^'^) - <I>f^5,,3,^ 



+ V2n - m - 1 max l^)?/:, n(e'^)l) 

which gives by ([80|) and ([HT]) . in conjunction with the facts that D{e^'^) / 
on [q,/3] and m(n) < n, relation ([83]) . □ 

5. Connection to positive qf on [-1,1] 

Remark 5.1. Let m,n G No,0 < m < 2n — l. Obviously, using the symmetry 
with respect to vr we have: 

y 

(85) / p{x)dij{x) = ^Xsp{xs) for p G \r2n~i-m, 

•' s=l 
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where — l<x„<...<xi<l and Ai, . . . , G if and only if 



n 



(86) / t{ip)da{ip) = Xs{t{(ps) + t{-ips)) for t G 7i„-i-m, 



where the ips = arccosx^ satisfy < ipi < < ■ ■ ■ < (pn < t^, and a is 
given by 

-0(1) — ip{cos 93) for < 99 < vr, 

ipi^cos If) — ^/^(l) for — TT < 93 < 0. 



Cj((/9) 



Note that, x = ^(z + i), 

n n 
s=l s=l 

By this equivalence and Theorem 13.11 we obtain immediately a complete 
description of positive qf of degree of exactness 2n — l — m on [—1, 1]. Indeed, 
Theorem 13.1b ) and d) will lead to characterizations given by the author in 
|16) . In fact in [16] we transformed the problem in the proof to the unit 
circle solved it there, that is, proved Theorem 13.1b ) and d) in the real case 
and transformed it back. It can be shown that the characterization by the 
recurrence relation, that is, by Theorem 13.1b ) yields, using the connection 
between the recurrence coefficients of OPUC's and the recurrence coefficients 
of polynomials orthogonal on [—1,1] (see [231 (13.1.7)] or [HJ Thm. 31.1]), 
the characterization of positive qf on [—1, 1] by the three term recurrence 
relation given in \17\ [19j by the author, see also [20j . 

When we consider Radau or Lobatto qf on [—1,1], see [201 Section 4] 
that is, that one or two nodes are the boundary points ±1 the correspond- 
ing equivalence is as follows: let us consider the Radau case, that is, the 
boundary point 1 is a node - analogously the other cases are obtained - then 

r " 

(87) / p{x)d^ix) = ^ XsP{xs) + Aop(+l) for p E Pan-m, 

s=l 

where — 1 < a;^ < • • • < a^i < 1 and Aq, Ai, . . . , An G M"^, if and only if 
7 A 

(88) / t{f)dai^) = y t(v3o) + Yl >^s{t{^s) + t{-(ps)) for t E T2n-m, 

— TT 

where the = arccos Xg satisfy ifo = 0<ifi<ip2<---<fn<'^ and 
where Aq, Ai, . . . , A„ E M"*". 
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